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ABSTRACT 

Using the Spitzer Infrared Spectrograph, we have performed mid-infrared spec- 
troscopy on the young binary brown dwarf 2MASS J04414489-h2301513 (15 AU) in 
the Taurus star-forming region. The spectrum exhibits excess continuum emission that 
likely arises from a circumstellar disk around the primary. Silicate emission is not de- 
tected in these data, indicating the presence of significant grain growth. This is one of 
the few brown dwarf disks at such a young age (~ 1 Myr) that has been found to lack 
silicate emission. To quantitatively constrain the properties of the disk, we have com- 
pared the spectral energy distribution of 2MASS J04414489-I-2301513 to the predictions 
of our vertical structure codes for irradiated accretion disks. Our models suggest that 
the remaining atmospheric grains of moderately depleted layers may have grown to a 
size of > 5 //m. In addition, our model fits indicate an outer radius of 0.2-0.3 AU for the 
disk. The small size of this circumprimary disk could be due to truncation by the sec- 
ondary. The absence of an outer disk containing a reservoir of small, primordial grains, 
combined with a weak turbulent mechanism, may be responsible for the advanced grain 
growth in this disk. 
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Introduction 



Grain growth, crystallization, and settling occur in primordial circumstellar disks around both 
stellar and substellar objects. These processes produce several observational signatures. The mid- 
infrared (IR) silicate emission becomes weaker as the grains in the upper layers grow and settle 
towards the midplane while crystallization modifies the shape of the band. IR continuum emission 
is reduced at longer wavelengths with grain growth and settling because the intercepted stellar flux 
decreases as the disk flattens. In the substellar domain, these hallmarks o f dust e volution have been 
reported previously for disks at ages t ^ I — 10 Myr (lApai et al.l 12004 l2005l : IScholz et al.l 120071 : 



Morrow et al.l 12008 : IPascucci et al.l l2009l : iRiad l2009l : iRiaz et al.l l2009l ) . For insta nce, disks around 
the young est brown dwarfs exhibit evidence of small, crystalline silicate grains (jApai et al.ll2005l : 
Riadl2009l ) whil e older brown dwa. r f disks have weaker s ilicate emission t hat indicates the presence 
of large grains (jScholz et al.l 120071 : iMorrow et al.l 120081 : iRiaz et al.l |2009| ) . The evolution of grains 
in the inner anr iuli of these disks, e.g., th e zone where the 10 //m silicate band emerges (within 
10~^ — 0.1 AU, iKessler-Silacci et al.l 120071 ) . appears to occur more rapidly than in disks around 
solar- mass stars. One possible explanation for this phenomena is that the low levels of turbulence 
for brown dwarf disks produce lower mass accr etion rates and reduce the replenishment of the inner 
disk with primordial grains from outer annuli (ISicilia-Aguilar et al.ll2007l ). Also, the collisional rate 
between large size grains would be small, preventing the fragmentation to small size grains. 

During a mid-I R spectroscopic survey of young brown dwarfs with the Spitzer Space Telescope 
([Werner et al.ll2004l ). we have identified a brown dwarf disk that shows evidence of advanced grain 
evolution at a relatively young age of ~ 1 Myr. This brown dwarf, 2MASS J04414489-h2301513 
(jLuhmanll2006l . henceforth 2M J04414489), lies in the Taurus star- forming region and was recently 
shown to harbor a companion with a projected separation of O'.'lOS 15 AU) and a mass of 5- 



10 Mjup (ITodorov et al 



2010l ). In this Letter, we use our mid-IR spectrum of this system and our 



models of settled, irradiated accretion disks to constrain the properties of the circumstellar disk 
that resides around the primary. 



Observations 



We observed 2M J04414489 with the Spitzer Infrared Spectrograph fIRS: iHouck et al.ll2004l ) 
on 2009 April 20 as a part of the Guaranteed Time Observations of the IRS instrument team. We 
collected data with both of the low-resolution IRS modules, Short-Low and Long-Low, which cover 
5.3-14 and 14-40 ^m, respectively, with a resolution of A/AA ~ 90. We used the basic calibrated 
data produced by version S18.7 of the Spitzer Science Center pipeline. Before extraction of the 
spectrum, bad and "rogue" pixels were fixed by interpolating adjacent, good pixels in the disper- 
sion direction. 2M J04414489 was then ex tracted using the optimal extraction tool in SMART 
(jHigdon et al.ll2004l : iLebouteiller et al.ll20ld ). We chose the default calibration, which is based on 
a relative spectral response function derived from three calibrator stars. The young low-mass star 
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2MASS J04414565+2301580 is 12'/4 from 2M J04414489 and fell within the slit during the LL obser- 
vations. We performed optimal extraction at the nominal nod position of 2M J04414489 to extract 
its spectrum, which should result in minimal contamination from 2MASS J04414565-I-2301580. In- 
deed, the resulting flux levels of the data at SL and LL match at 14 ^m without applying any 
scaling factors. The IRS spectrum includes the flux of both 2M J04414565 A and B, but it should 
be dominated by the former, as explained in Section 13.11 



3. Analysis 
3.1. Disk Emission 

We show the spectral energy distribution (SED) of 2 M J04414489 in Figure [H which was con- 
struc ted with our IRS data , 0.8-2.5 //m spectroscopy from lLuhmanI (|2006l ). and 3-8 ^m photometry 
from lLuhman et al.l (j201Cl ). We have included in Figure [1] the SED for a diskle ss young brown dwar f 
near the same spectral type as 2M J04414489 (2MASSW J1139511-315921, iMorrow et ahlbood '). 
Relative to the stellar photosphere of the latter object, the SED of 2M J04414489 exhibits significant 
excess emission that indicates the presence of a dusty, optically thick disk, which is consistent with 
earlier results based on the mid-IR photometry (jLuhman et al.ll201Cll ). The disk is probably located 
around the primary because a disk around the secondary would be too faint to acco unt for the ob- 
serve d excess emission based on the flux ratio at near-IR wavelengths (AX^ 



1.54, 



Todorov et al 

2ninl ). 



20 id ) and the typical near- to mid-IR colors of brown dwarfs with disks ( Luhman et al.ll20ld ). A 
circumbinary disk is not plausible since its inner hole would be too large for the disk to produce 
significant mid-IR emission. Finally, the mid-IR excess emission cannot be due to the photosphere 
of the secondary since its IR colors should be similar to those of the primary's photosphere. For 
instance, given the spectral type of M8.5 for the p rimary and an expected type of M9.5-L0 for 

20ld). their ph otospheres should have colors near 



2006 



Todorov et al 



the secondary (jLuhmanI 
Kg — [4.5] = 0.71 and 0.83, respectively (iLuhman et al.ll20ld ). Combining those colors with the 
fiux ratio at Kg suggests that the secondary's photosphere produces an excess of only -^0.02 mag 
in Kg — [4.5] for the system relative to the photospheric colors of the primary. 



3.2. Model Parameters 



The optical spectrum of 2M J04414489 exhibits Ha emission ([Luhmanl 120061 ) . indicating the 
presence of ongoing mass accretion. Therefore, we have modeled the mid-IR excess emission from 
2 M J04414489 in t e rms of emis s ion fr om a settled, irradiated accretio n disk followin g the procedures 
of b'Alessio et al.1 (|l998l . bood . bood l adapted for substellar objects (jAdamelboim . For the steflar 
photo sphere, we adopt an effective tempe rature of 2555 K, which corresponds to its M8.5 spectral 



type (jLuhman et al 



2003 



Luhmanl 120061 ). We also adopt a luminosity of 0.004 Lf^ based on its 



J magnitude and a mass of 0.025 Mq based on theoretical evolutionary models (IChabrier et al 
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2000) • In our models, we solve the vertical structure of the disk while assuming a constant turbulent- 
viscosity parameter a and a uniform mass accretion rate, M, through the disk. For most of the 
brown dwarfs in Taurus, assuming a > 10~^ is suitable ( Adamd I2OIOI ). although we allow a to 
range from 10~^ to 0.01 for our models of 2M J04414489. In each annulus, the mass surface density 
depend on the ratio M /a. We compute models with M = 10~^^, 10~^\ and 10~^^ M.(^ yr~ ^, which 



are typical values ni e asured for young substellar objects (jNatta et al.ll2004l : iMuzerolle et al 



Mohanty et al 



2005 



Herczeg &: Hillenbrand 



2008 



Herczeg et al.ll2009l ^. 



2005 



The adopted dust model is composed of seg regated spheres of a stro nomical silicates and 
graph ite with abundances and optical constant s fromlDraine fc Led (| 19841 ) and lWeingartner &: Draine 
(|200ll ) and a size distribution of n{a) ^ a~P ([Mathis et al.l 119771 ) . where p = 3.5. To account for 
grain growth, we have varied the maximum grain size in the upper layers and the inner wall, adopt- 
ing values of Umax = 0.25, 1, 5, 10, 100 and 1000 /im. Meanwhile, we assume Umax = 1 mm for 
the dust population near the midplane. As grains grow and settle towar ds the midplane, the at- 
mosphere is depleted of small grains. The model of lD'Alessio et al.l ()2006l ) quantifies this depletion 
through the parameter e, which is the dust-to-gas ma ss ratio of the atmo sphere relative to the 
standard ratio (with = 0.004 and Cgraph = 0.0025, iDraine Sz Ledll984l ). We have considered 
e = 0.001, 0.01, and 0.1 in our models. 

We also calculate the str ucture and SEP of t he inner wall , which is heated by radiation 
from the stellar photosphere (jD'Alessio et al.l l2005l ) . The wall is located at radius Rin , which 
can be the dust destruction radius (assuming a sublimation temperature of ~ 1400 K) or the 
magnetospheric radius, whichever is larger. We explore maximum dust temperatures of ~ 1400 K 
to 1000 K (in steps of 100 K), corresponding to an inner dust gap with Rjn, ^ 3 to 6 R^. . The 



outer disk radius Rd may b e set by interaction with the secondary (jPapaloizou fc Pringld 119771 : 
Artvmowicz fc Lubowl Il994l ). Assuming a mass ratio of the s econdary to the primary of ~ 0.3, 
a projected semimajor axis of ^ 15 AU (jTodorov et al.ll2010l ). and a mean Reynolds number of 



truncation ^ 



10^-10^ f see I Artvmowicz fc Lubowlll994l ). the truncation radius is predicted to be 3 < iZ, 
7 AU for eccentricities between e = and e = 0.5 — 0.6 (a large circumprimary radius corresponds 
to a small eccentricity). Therefore, we fixed Rd to an initial value of 5 AU; for promising models, 
we have selected values of Rd that range from 0.1-5 AU. Finally, we adopt inclination angles of 
i = 10, 20, and 60°. The disk SED library consists of some 700 models. 



3.3. Model Fits 



We proceed to identify the combinations of disk parameters that can reproduce the SED of 
2M J0441 4489. We perform a re duced-y^ goodness-of-fit test between each model and the IRS 
data (e.g. Ivan Boekel et al.l l2005l ) . where ~ 1 indicates a good match. We find that models 
with z = 20°, 10"^^ < M < 10-1° M© yr'^, and a wah temperature of ~ 1000 K match the 
SED at A < 10 nm better than other values of these parameters. The radius of this inner wall is 
Rin ~ 6 R* ~ 0.009 AU. We allow the wall height z^aii to vary from 3.4 Hp to 4.1 Hp, where Hp is 
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the gas scale height at The models with these values of i, M, and Rin that produce Xr < 4 
are 21. An important constrain fo r our model SEP selecti o n is t he value of the mass accretion rate. 
From the M — M^, relationship of lHerczeg &: Hillenbrand! (|2008l ). the mass accretion rate onto this 
low mass brown dwarf would be logM ~ —10.7 it 0.6, thus models with 10"^'' yr~^ are in the 
higher end of the expected mass ac cretion rate. O n the other hand, the equivalent width of the 
Ha emission hue of 2M J04414489 (ILuhmanl 20061 ') is 100 A, which corresponds to a 10% width 



of Ha of ~ 220 km s"^ (Muzerolle et al 



20051 ). and to a mass accretion rate of < 10 Mq yr 



(jNatta et al.ll2004l ). Hence we discard the high mass accretion models, and select the remaining 
four models as our best-fit models, which are listed in Table [TJ 

The SED predicted by model A, which has the lowest value of Xr; is compared to the observed 
SED in Figure [H The variations in the disk parameters among the best-fit models in Table [T] are 
qualitatively explained by the fact that the mid-IR emission is weaker for smaller e and M and 
larger a^ax and a, causing a reduced optical depth. For instance, the grain growth of model C and 
the depletion of model D must be balanced by increasing the dust mass content in order for the 
models to produce sufficient emission to match the data. Meanwhile, models with a^ax = 1 miii 
and most of the models with 10~^^ < M < 10^^^ Mq yr^^ produced emission that was too weak 
to reproduce the data (x^ ^ 4). 

We compare the silicate emission predicted by our best-fit models to the IRS data for 2M J04414489 
in Figure [21 We normalize the IRS spectrum and the models with a continuum interpolated from a 
third-order polynomial constru cted by fluxes between ~ 5.2 and ~ 7.9 /xm and between ~ 13 and 



16 ^m. (e.g. iFurlan et al.ll2006l ). The strength of the silicate feature relative to the continuum is 



obtained by subtracting the continuum to the feature, int egrating this residu al flux and normalize 
it to the integrated continuum flux between 8 and 13 /xm ([Furlan et al.ll2006l ). The strength of the 
observed feature is 1.03, whereas the strengths of our models features range from 1.024 to 1.032. 
Thus, the observed and modeled silicate features are similar, but also, they are extremely weak. 
The absence of detectable silicate emission in the data indicates that grains in the disk atmosphere 
have grown to sizes of 5 < dmax < 1000 fiui, which are larger than the typical ISM grains. 

The dust mass es of our final set of disk models are listed in Table [H using (^^std = 0.0065 
(jPraine Ledll984l ). The masses of the models range from ~ 3 to ~ 5 Mmooti, which are more 



than an order of magnitude lower t han the values es timated from millimeter measurements of 
other brown dwarf disks in Taurus (jScholz et al.l 120061 ). This result is not surprisi ng given that 



close binary stars exhibit low disk masses compared to single stars and wide binaries (jJensen et al 
19941 ^. 



3.4. Size of the Mid-IR Emission Region 



The mid-IR emission from an isolated disk irradiated by a M8.5 brown d warf typically forms 
within Rmir. ^ 5 AU. The exact size of this region depends the size of the grains (iKessler-Silacci et al, 
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20071 ). but also on the dust mass surface density through the variation of e, M or a. We find that 
decreasing a or increasing M or e result in a wider emitting region since those variations increase 
the irradiation surface and the impinging irradiation flux at any given radius. On the other hand, 
decreasing amax will increase the optical depth to the incident radiation at the upper layers, which 
in turn extends vertically outward the irradiation surface, allowing more distant annuli to contribute 
to the emission. 

If the outer radius of the disk is smaller than the nominal radius of the mid-IR emitting region, 
then the mid-IR SED is sensitive to the value of the outer radius. For instance, the radii for all of the 
models in Tableware constrained to be 0.2 < Rd < 0.3 AU by the mid-IR SED of 2M J04414489, 
which is smaller than the expected radius given by the binary tidal interaction {Rtruncation ~ 3-7 
AU). However, if the eccentricity of the binary is high enough, the periastron would be close to 
1 AU if a semi-major axis of 15 AU is still assumed. The dependence of the predicted SED on 
the outer radius is illustrated in Figure [3] for model A. Although the disk radii for our models of 
2M J04414489 are roughly similar to that derived for the circumprimar y disk of SR 20 (0.39 AU, 



for a G7 primary and a binary projected separation of 5.3 — 9.9 AU, iMcClure et al.M2008l ). the 
former has a redder SED (i.e., less affected by truncation) because the mid-IR emitting region of a 
brown dwarf is much smaller than that of a star. 



Discussion 



We find that the dust in the disk around 2M J04414489 is relatively evolved for the age of the 
system (~ 1 Myr). Based on the slope of its SED and the absence of silicate emission, most of the 
initial atmospheric grains have settled towards the midplane (0.01 < e < 0.1) and the remaining 
grains may have experienced significant growth (5 < dmax < 1000 /xm), representing one o f ^ 2 
brown dwarf disks in Taurus that show evidence of accelerated dust evolution (see KPNO iRiaz 
20091 ). Th is degree of dust evolution has been observed primarily for brown dwarfs at ages of 
5-10 Myr (jScholz et al.ll2007l : iMorrow et ahlboosl l. 



Some mechanisms to explain grain grow th and settling in di sks around either single or multiple 
late-type M-dwarfs have been put forward. iMeeus et al.l (120031 ) suggested that binary interaction 
may accelerate the dust evolut i on th rough an enhanced vertical stirring of grains from the mid- 
plane. However, iPascucci et al.l (|2008l ) and Furlan et al. (2011, in preparation) found no significant 
differences between the silicate feature for early type single and binary star systems in Taurus, re- 
vealing that binary interaction may not affect the inner disk region, althoug h the outermost region 
is indeed influenced. The interac tion truncates the disk to a given radius ([Papaloizou &: Pringle 
19771 : lArtvmowicz &: Lubowlll994l ). decreasing the emission from a certain wavelength range. For 



^For the other objects reported by Riaz 1 20091 ) . siHcate emission is detected at weak levels in the same spectra 
that were re-reduced (Furlan et al. 2011, in preparation). 
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instance, the submillimeter and millimeter emissi on from the circumprima ry disk is diminished for 
projected binary separations of ~ 50 — 100 AU (jJensen et al.lll994l . Il996l ). while for smaller sep- 
arations (< 50 A U), the mid-IR emiss ion of the outer disk decreases with respect to that of the 
innermost region (jMcClure et al.ll20ld ). The mid-IR emission of J04414489 is the result of this 
truncation, whereas the fluxes near ~ 10 /im are unlikely influenced by the binary interaction. In- 
stead, the inner region may be probably affected by the turbulence as most of the eddie dissipation 
occurs within those annuli (< 0.1 AU), which represents ~ 30 — 50% of the total radial extension 
of the disk. 

Grain growth and settling should be efficient in regions close to the star, where tempera- 



tures and densities are larger, and the turbulent move ments are higher (jBeckwith et al.l I200C 



Dullemond fc Dominik 



2004 



Nomura &: Nakagawal l2006l . references therein). For this weakly tur- 



bulent disk (a = 10~^), the coagulation rate may be slowed down, but also the relative velocities 
of the grains, thus preventing their fragmentation. The smaller velocity gradients can work to pre- 
serve the large grains once formed. A radial mixing mechanism may exist that moves evolved grains 
outward and replenishes the inner disk with primordial dust (e.g. turbulent diffusion, conve ction. 



gravitational torques, among others, see [Gail 



2001 



Bockelee-Morvan et al 



2002 



Bosdl2004l ). The 



primordial dust reservoir is limited by the outer radius of a disk. If a disk is outwardly truncated to 
a small enough radius (e.g., Rd < Rmir), the outermost region will contain evolved grains, which 
return to the innter region to be processed f urther. The characteristic tim escale for this radial 
mixing is approximately the viscous timescale (jBockelee-Morvan et al.ll2002l ). which for an a-disk, 
tRM = Ro^K{Ro)/'icicl{Tc). Here, Cs{Tc) is the sound speed at the midplane and Q,k the Keplerian 
angular speed, both evaluated at Rq (the initial outer radius of the disk). We evaluate Irm at the 
outer radii listed in Table [T] assuming a central temperature at the midplane of ~ 60 — 80 K. The 
mixing timescale derived is ~ 1 — 2 Myr, which is comparable to the age of 2M J04414489. The 
radial mixing may be at an earliest stage, but its influence over the dust processing cannot be 
neglected. Therefore, low levels of turbulence combined with a highly truncated disk could explain 
the deficit of primordial dust that we observe in the disk of 2M J04414489. 



5. Conclusions 

We have presented a mid-IR spectrum of the binary brown dwarf 2M J04414489 in the Taurus 
star-forming region. The spectrum exhibits excess emission that we attribute to a circumprimary 
disk, but silicate emission at 10 fim is not present, indicating that significant grain growth has 
occurred, leading to a disk on which small grains no longer exist (5 < cimax < 1000 /um). This is 
one of ~ 2 brown dwarf disks at such a young age (~ 1 Myr) that have been found to lack silicate 
emission. Our models of the SED of 2M J04414489 suggest that the outer radius of its disk is rather 
small {Rd ~ 0.2-0.3 AU), which may reflect truncation by the binary companion. The absence of 
an outer disk with a reservoir of primordial grains, combined with weak turbulence, could explain 
the advanced grain growth in the inner disk of 2M J04414489. 
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Fig. 1. — SED of the young brown dwarf 2M J04414489 compared to a model of its circumstellar disk 
(dotted and dashed lines, mo del A from Tabled]). We have adopted the SED of 2MASSW J1139511- 
315921 (iMorrow et al.ll2008l ) to represent the stellar photosphere and have scaled it to the 1-2 /um 
flux of 2M J04414489. For A < 2.5 /im, the stellar photosphere overlaps almost precisely with the 
photosphere of 2M J04414489. 
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Fig. 2. — Spectrum near 10 /xm for 2M J04414489 normalized to the continuum surrounding the 
wavelengths of the silicate feature. For comparison, we show the normalized silicate emission 
produced by the models in Table [TJ The order of the models from strongest to weakest emission 
peak is B, A, C, and D. 
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Fig. 3. — SED of the young brown dwarf 2M J04414489 compared to model A (Rd = 0.2 AU) and 
two models that have the same parameters as model I except for different values of Rd- 
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Table 1. Models^ 



Model^ 


€ 


dmax 




Rd 


Mdust 




tRM 








(Hp) 


(AU) 






(10^ yrs) 


A 


0.1 


10 


3.7 


0.2 


3.2 


3.7 


1.2 


B 


0.1 


5 


3.6 


0.2 


3.2 


3.8 


1.2 


C 


0.1 


100 


3.9 


0.3 


5.2 


3.9 


1.8 


D 


0.01 


5 


3.8 


0.3 


5.7 


3.9 


1.9 



^The listed parameters are: the dust-depletion paramenter e, 
the maximum grain radius of the dust population in the disk at- 
mosphere ajnaxj the height of the wall z^aU, the outer disk radius 
Rd, the total dust content of the disk Mdust, the reduced-x^ value 
of the model, and the characteristic timescale for the radial mixing, 
evaluated at Rd, t^M- 

= 10-^ M = 10"^^ Mq yr-\ i = 20°, Rin = 6 for aU 
models. 



